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Abstract. We present simulations of the dynamic critical state for a 2D superconductor with strong pinning 
centers, corresponding to a matching field twice the applied magnetic field. Fhe sharp crossover between 
the plastic regime, at low current density and temperature, and the fluid flow regime for flux motion is 
characterized by the activation energy for flux motion and the transverse diffusion of the vortices trajectory. 

PACS. 64.60.Ht Dynamic critical phenomena - 74.60.Ge Flux pinning, flux creep and flux-line lattice 
dynamics 

1 Introduction non zero temperature. However, numerical simulations [[j} 

0,0,0], experiments |^,^| and theory have shown that 
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^3 Flux pinning in type II superconductors is a well docu- a driven flux lattice can be a rather complex object, re- 
Ch ! 

O mented example for the general problem of an elastic pe- mar kably different from the simple pictures used in the 
! riodic medium interacting with disorder. The study of the interpretation of the experiments. In particular, the pin- 
dynamics of the pinned flux line lattice, under the action n i ng f a two-dimensional flux lattice by strong pinning 
of an external current or by the virtue of a time varying centers lends itself easily to numerical simulations and has 
external magnetic field, is one way to apprehend the prob- been the subject of several studies. Such a system is also 
lem of flux line pinning. Well known examples are trans- a goo d approximation for highly anisotropic superconduc- 
port I-V characteristics and magnetization relaxation ex- tors irradiated by heavy ions, such as the high Tc oxides, 
periments. Such experiments are usually interpreted using wn i c h are now well documented. In this paper, we present 
simple models, assuming a uniform critical current den- numerical simulations of this strong pinning situation, at 
sity and a single, current dependent, activation energy at 
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non zero temperature and for conditions similar to those 
encountered in experimental magnetization experiments. 

2 Simulation 

We simulate a slab submitted on one side to an external 
magnetic field, applied parallel to its boundary, while, on 
the opposite side, the external magnetic field is assumed to 
be zero (vortices are removed as they cross this boundary). 
Periodic conditions are used, so that the effect of the finite 
dimensions is limited. Flux lines are assumed rigid rods. 
This correctly modelizes a layered superconductor with 
decoupled layers in the critical state, for which the local 
magnetic induction decreases from the applied field value 
on the external sides to a lower value at the center of the 
sample (Fig. Q). 

Strong pinning centers are randomly distributed in the 
sample. The density of the pinning sites is n = B$/<I>q, 
with the flux quantum and B$ the 'matching field' for 
which an equilibrium flux line lattice shows the same den- 
sity of flux lines. A region of width A from the side where 
vortices are introduced is kept free of pinning sites, in or- 
der to allow for the initial formation of a regular hexagonal 
lattice in this region. The pinning sites are assumed nor- 
mal cylinders parallel to the applied field, with radius cq. 
The situation where, at low temperature, the vortex core 
radius £(T) is smaller than cq and pinning is due essen- 
tially to the reduction of the core energy when the line 
sits on the pin is considered. The force exerted by a pin 
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at a distance r to the line is given by: 

{v£ o r/r £ for r < r , 
(1) 
for r > ro . 

where £o = ($o/47rA) 2 is the line energy, r = c + £/2 
and a < 1. 

The force per unit length exerted by a vortex at a 
distance r to the line is: 

Mr) = K/8tt 2 A 3 ) Ki (r/A) (2) 

where K\ is a Bessel function. This is a good approxima- 
tion strictly only in the case of vortex lines (rods) and 
for 2D vortices a logarithmic interaction should be used. 
In the present case, the more rapid decrease of the Bessel 
function allows us to cut the interaction between vortices 
at a distance 5A and save computation time. 

The external magnetic field B$ is simulated by an ex- 
tra force fso acting on each vortex, perpendicular to the 
external side of the sample. The force acting on a vortex 
at a distance x from the boundary is the one imposed by 
a semi-infinite vortex lattice at a distance ao -I- x, where 

1/2 

ao = {$o/B) ' is the flux lattice spacing at the equilib- 
rium. 

The finite temperature is simulated in a way similar to 
the one described in ref. jj| by adding a stochastic velocity 
to the particles: 

v T - (2 k B T r/3Z\ 2 ) 1/2 ( 1X I 0(At/r - q) (3) 

where 1 jr is the average frequency for the thermal pertur- 
bation, 7 Xj j, are random numbers from a Gaussian distri- 
bution of width 1, q is a random number from a uniform 
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distribution between and 1, and O is the Heaviside func- 
tion. We have used a single set of parameters r and At 
for all temperatures with r> /li. The accuracy for the 
temperature simulated in this way was checked by com- 
puting the life time of a single vortex in a potential well 
similar to the ones defined above. The life time as a func- 
tion of the well depth was found to match accurately the 
one given by the Arrhenius-Kramers relation for an over- 
damped particle, in the range of temperature studied here. 
In this overdamped regime, results are independent of the 
viscosity when expressed in the time unit t/rj 1 the escape 
frequency being inversely proportional to r\. 

Finally, the displacement Av\ of vortex i during time 
step At is computed according to the diffusive equation: 

T) AvJAt = 

f w(^) + Y { p( r v) + f B (^/ao) + V v T (4) 

j¥» 3 

where the first sum is over all vortices at a distance less 
than 5A and the second one is over all pins (77 is the vis- 
cosity per unit length). The time step is chosen so that 
the vortex displacement during one step is always much 
smaller than the pinning characteristic length, ro- At non 
zero temperature, vortices tend to creep through the sam- 
ple from the boundary submitted to the external field, 
where vortices are fed whenever the magnetic pressure is 
strong enough to allow for their penetration, to the one 
with zero applied field. As a consequence, a steady ther- 
mally activated flux flow is set after some time, charac- 
terized by a linear flux density profile and a steady vortex 
current density, T , that flows through the sample. 
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3 Results 

In the rest, we use the following material parameters, ap- 
propriate for a high temperature superconductor: A = 
1400 A, £ = 18 A, c = 35 A and a = 0.1. The sam- 
ples held typically 5000 pinning sites and 1000 vortices in 
the steady state. The matching field, B$ — 5000 G and 
the external field, Bq = 2500 G, were the same for all the 
simulations. 

In a way similar to what is done in the analysis of real 
relaxation experiments, we use a one-dimensional model 
with a single activation energy to analyze the vortex cur- 
rent density JlC||. According to this widely used model, 
vortices thermally hop over barriers disposed along a 1- 
D axis with height U(j) (where j — (47r) _1 9i? / 'dx is 
the current density) and the vortex current density is 
T = j B 0o V 1 ex P(—U(j)/kBT). This model clearly 
misses some of the fundamental features of the real ex- 
periments. As pointed in ref. Jl0| , it is valid only in the 
case of almost uniform vortex and current densities, when 
physical quantities may be replaced by their average over 
several vortex spacing. In particular, it assumes that the 
screening current flows in our case along straight lines per- 
pendicular to the vortex flow. As can be seen in Fig. ||A, 
this is not so and the current tends to flow along some 
curved, branched paths that reveal the highly inhomoge- 
neous stresses in the pinned flux lattice. There has been 
attempts to account for the spatial heterogeneity of the 
real systems, using a distribution of activation energies [pT[ 
|l2| |. However, the method assumes parallel, independent 
relaxation channels, which is still a very crude assump- 
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tion. As for real magnetization experiments analysis, we 
obtain the average screening current density from the to- 
tal magnetization of the sample. In order to investigate 
different flux gradient profiles, different sample geometries 
were used with the length of the cylinder between 20A and 
I00A and diameters between 10A and 300A. The experi- 
ments probed different T and j values, which are located 
on a smooth (T, j) trajectory shown in Fig.^j. This is sim- 
ilar to relaxation rate measurements in superconductors, 
for which both the temperature and the screening current 
are varied - the latter parameter being implicitly deter- 
mined by the temperature and the time window explored 
in these experiments. 

The activation energy obtained from the simulations 
is shown in Fig. |^, where energy is normalized to the sin- 
gle vortex pinning site energy, Uq — aeo^o/^, and crit- 
ical current density to the single vortex critical current, 
jo = (7£o/£0o- It is shown as a 3D plot, as the computa- 
tion of the energy with one of these two parameters fixed 
would imply prohibitive computation time as temperature 
or current decreases. The same limitation is also encoun- 
tered in real relaxation experiments, although in this case 
the temperature variations are usually disregarded and 
the data plotted as U(j). The reason for this is that Ma- 
ley's procedure Jl3| allows to show that, in most cases, the 
activation energy obtained from relaxation experiments is 
temperature-independent. We notice, however, that this 
is no longer true at low temperature, which has been 
attributed to the occurrence of the quantum tunneling 
regime Q. 
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A clear upturn can be seen in Fig. |^, where the activa- 
tion energy increases strongly with decreasing current be- 
low j/jo — 0.05, ksT/A ~ 0.25. This upturn corresponds 
to the crossover between a plastic to a fluid flow regime 
for flux motion. This can be seen, first, from a direct ob- 
servation of the vortices trajectories: in the plastic regime, 
there exists a time scale at which some vortices move over 
distances larger than do, while others remain immobile 
on that scale and constitute some kind of 'pinned islands' 
(Rcf. Q). It is important to stress that these islands de- 
pin at a larger time scale for non zero temperature, so 
that only a dynamical definition of these domains may 
be given in our case. No such characteristic time may be 
found in the fluid flow regime where, depending on the 
time scale considered, nearly all or none of the vortices 
move over this distance (Fig. §A,B). A more quantita- 
tive argument in favor of two clearly distinct regimes is 
given by the examination of the diffusion of a vortex tra- 
jectory, transverse to the average flux flow (Fig. |^ and 
^|). In the fluid flow regime, the transverse displacement 
on the scale of the equilibrium flux lattice parameter, 
S = O.Q 1 ([y(0)~ y(ao)] 2 ) 1 / 2 where the average is performed 
along the vortex trajectory, is significantly smaller than 1, 
whereas it increases with decreasing temperature and in- 
creasing current density up to a value ~ 0.6 — 0.65 in the 
plastic regime (Fig. ^). This limiting value may be under- 
stood by considering that the most disordered trajectory 
allowed in the flux lattice is the one where a vortex hops 
randomly from one interstitial site to a neighbor site. A 
consequence of disorder in the plastic regime is that the 
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examination of the flux pattern hardly reveals the orien- 
tation of the average flux flow (Fig. ||C, making abstrac- 
tion of the region close to the zero field boundary of the 
sample, where side effect is visible), as compared to the 
fluid flow regime (Fig. ||B). As can be seen in Fig. ||, the 
increase of the transverse displacement along the (T, j) 
trajectory shown in Fig. || is rather abrupt, suggesting 
that the trajectory might cross a transition line rather 
than a crossover line. If this is so, an order parameter 
may be elaborated from 5, after subtraction of the high 
temperature baseline and adequate normalization, which 
is unity in the ordered regime and tends towards zero in 
the opposite limit, thus delimiting two distinct, nonequi- 
librium, phases. The static topological orders of the two 
phases which would be defined in this way do not show 
significant differences: both lack long range order, either 
crystal-like or smectic-like. To that respect, they both look 
like a liquid, as demonstrated by the inspection of the den- 
sity autocorrelation function (Fig. ^|). It is difficult, how- 
ever, to unambiguously put into evidence the existence 
of a second-order noncquilibrium phase transition in our 
case. 

Up to now, we have reported only upon a single point 
of the crossover line in the (T,j) plane. A close examina- 
tion of Fig. H provides some indication about the general 
behavior of this line. Indeed, it can be seen that both sam- 
ples B and C do not exhibit a strictly homogeneous regime 
and that the region close to their bottom boundary tends 
to be in the fluid flow regime. This is because, due to edge 
effects, the average screening current density tends to be 
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larger in this region. As a consequence, we may now pro- 
pose a schematic phase diagram for flux motion, where 
we have added the disorder strength as a third dimension. 
These observations are not easily compared to the theoret- 
ical expectations for driven lattices ([|| and refs therein). 
First, models investigate the case of a weak random poten- 
tial, which is not adequate here, unless the matching field 
greatly exceeds the applied magnetic field and a <C 1. 
Then, the two dimensional, non zero temperature case 
for the random potential is far less documented than the 
three dimensional oneQ. We may simply notice that the 
measured activation energy at the crossover is close to 
the plastic energy barrier (i.e. the barrier experienced by 
a moving interstitial defect) U p i ~ £o/2\/37r ~ 0.4 Uq 
(Fig. ||) Jl(| and that it is likely that this condition holds 
all along the crossover line. Indeed, energy barriers larger 
than U p i cannot be allowed by decoupled vortex flow chan- 
nels and larger barriers can be sustained only by the cou- 
pling of these channels in the fluid flow regime. Finally, we 
would like to compare with real relaxation experiments. 
The main difficulty in realizing the situation described by 
our simulations comes from the presence, in real crystals, 
of naturally grown point defects. Indeed, as-grown crystals 
invariably exhibit at low temperature high critical cur- 
rent. In the highly anisotropic Bi22i2 material, screening 
current as large as 10 10 AmT 2 is observed at low temper- 
ature. The introduction of strong pinning centers, such as 
columnar defects, does however induce a sizeable increase 
of flux line pinning at low temperature (almost one order 
of magnitude for the critical current density at 4.2 K in 
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Ref.JlT]], at a dose = 20 kG). From this, it may be con- 
cluded that natural defects are weak or dilute enough, so 
as to be overcome by the introduction of extra columnar 
defects, and that irradiated crystals might provide a situa- 
tion close to the one described here. An expected signature 
for the crossover is a weaker temperature dependence of 
the (relaxed) screening current in the fluid flow regime - 
as compared to the plastic one - due to the increase of the 
energy barrier in this regime. A sharp change is indeed 
observed in the temperature dependence of the relaxed 
screening current, for both as-grown and irradiated sam- 
ples, although this has been interpreted as the signature 
of the crossover between small and large bundles regime 
in the first case, and depinning from the defects in the 
second onejlTj]. We believe that a decisive test should be 
the measurement of the energy barriers at this point and 
the comparison with the one for the plastic mechanism. 

4 Conclusion 

We have studied the crossover between the fluid flow and 
the plastic regimes of the critical state dynamics for a two 
dimensional flux lattice, with strong pinning sites corre- 
sponding to a matching field twice the applied field. Fol- 
lowing a smooth trajectory in the (T, j) diagram, we find 
that the transverse displacement of the vortex trajectory 
rises abruptly as this trajectory crosses the crossover line. 
There, the apparent activation energy for flux motion is 
close to the plastic energy barrier and increases strongly 
in the fluid flow regime. We point out strong similarities 
between these results and the sharp crossover observed 
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in the temperature dependence of the screening current 
density for as-grown and irradiated Bi22i2- 
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ted line indicates the plastic barrier, U p i. Labels 'B' and 'C 
refer to the samples in Fig.|[ 

Fig. 3. Sample A: static 'critical state' at T = 0. Lines are 
colinear to the local current direction and their length indicate 
the current intensity (non linear scale). Samples B and C: same 
notation as in Fig. ^; lines are vortices trajectories in the steady 
state. 

Fig. 4. Displacement transverse to the average velocity, along 
vortices trajectories, j/(i)(the average vortex flow is along the 



Fig. 5. The transverse displacement on the scale of the equi- 
librium lattice parameter. Lines are guides to the eye. 

Fig. 6. Left: in gray scale, density autocorrelation function 
of the lattice snapshots. Right: after radial integration. Top 
and bottom are respectively sample B and C in Fig. [|, both 
showing liquid order (The exponential decay length, as defined 
in Ref. [jB"| , is respectively 0.8 ao and 0.5 ao). 

Fig. 7. Schematic phase diagram. The arrow indicates the 
trajectory followed in Fig. ^[ with indication of samples B and 
C in Fig. | and | . 
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